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ABSTRACT 



Harwit, Protheroe, and Biermann (1999) recently proposed that Bose- 
Einstein photon bunching might significantly affect the interpretation of 
Cerenkov counts of TeV gamma photons. Here, we show that a combination 
of two recent results of Aharonian et al. (2000) and Aharonian et al. (2001) 
permits us to set new, more stringent upper limits of <^ 10% on the fractional 
amount of photon bunching in the 7-10 TeV radiation from Markarian 501. 
Potential bunching at even higher energies should nevertheless continue to be 
investigated for this and other TeV sources, since a clear understanding of TeV 
energy spectra is required to unambiguously determine the spectral energy 
density of the mid-infrared extragalactic background. 



Subject headings: gamma rays: observations; BL Lacertae objects: individual 
(Markarian 501); radiation mechanisms: general; infrared: diffuse radiation 
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1. Introduction 

Harwit, Protheroe and Biermann (1999, hereafter Paper I) have proposed that 
some of the highest energy TeV Cerenkov gamma-showers attributed to single photons 
arriving from distant blazars might actually be produced by bunched photons at lower 
energies. Lower energy photons can penetrate greater distances through extragalactic space 
without suffering annihilation through collision with photons from the intergalactic infrared 
background and producing electron-positron pairs. The bunching proposed in Paper I was 
just the ordinary Bose-Einstein bunching well known at lower energies. Levinson (2000, 
2001) has recently re-analyzed this possibility, basing his study largely on a laser mechanism 
which, however, obeys very different photon statistics. This letter provides a more detailed 
analysis of photon statistics, and goes on to discuss the specific properties of the TeV 
radiation from Mrk 501. Section 2 of this letter provides a quantitative treatment of photon 
fluctuations and correlations from incoherent sources. Section 3 sets new upper limits on 
the Bose-Einstein bunching of TeV photons from Mrk 501. Section 4 touches on photon 
fluctuations in coherent processes like the laser mechanisms investigated by Levinson (2000, 
2001), and shows how these differ with specific reference to TeV gamma photons. Section 
5 stresses the importance of more stringent TeV observations in order to better define the 
infrared extragalactic background. A final section summarizes the main conclusions. 

2. Quantitative Treatment 

The photon fluctuations in a beam of radiation can be derived in a number of different 
ways (Hanbury Brown & Twiss, 1957, 1958a, 1958b, 1958c; Purcell, 1956; Wolf, 1957; 
Mandel, 1958; Harwit, 1960). Let us assume first that the radiation is blackbody at some 



-4- 



temperature T. Then, the Einstein-Fowler equation tells us that 



((AE) 2 ) = kT 2 ^ = hv(E) 



1 + 



1) 



exp(hu/kT) 

where the brackets () denote ensemble averages, T is the source temperature, E is the 
energy of the individual photons and v is their spectral frequency. From this we can 
immediately obtain the fluctuation ((AiV) 2 )in photon number N in a given volume and its 
relationship to the mean value (N) 

1 



((ANY) = (N) 



1 + 



(N) 



1 + (N)/g 



exp(hiy/kT) — 1. 

where g is the number of phase cells in the volume occupied by the photons. 



(2) 



g = 2AtcQvoAvc~ 3 (3) 

and A is the aperture of the telescope, r is the time interval over which a single measurement 
in the ensemble of measurements of N is made, Q is the solid angle subtended by the 
source, u is the central frequency and Av the radiation bandwidth. The numerical factor 2 
accounts for the two possible states of polarization. 

If the radiation incident on the light-gathering telescope is coherent, in the sense that 
the radiating source just fills the telescope's diffraction angle at frequency u , the transverse 
dimension of the incoming beam is simply the cross sectional area of a single phase cell 
and, for fully polarized photons, the number of phase cells along the direction of the beam 
is g ~ tAv. Though we have derived this relationship for blackbody radiation, Mandel, 
Sudarshan and Wolf (1964) showed that it applies to any quasi-monochromatic coherent 
beam of radiation, regardless of whether or not it originates in a statistically stationary 
thermal source, as long as the wave amplitudes in the beam are represented by Gaussian 
random variables (Mandel, 1963). 

Equations (1) and (2) specifically, do not apply to lasers or masers operating well above 
their thresholds; this is further discussed in section 4, below. 
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In practice, we do not directly measure numbers of photons, but rather instrumental 
counts produced by the photon flux. If the counting efficiency of the instrumentation is e, 
the fluctuation in the number of photoelectrons, N t , becomes 



where the values of e and g characterizing the instrumentation and the incoming stream 
of radiation may be experimentally determined. This permits us to calculate the true 
fluctuations in the incoming radiation beam {(AN) 2 ) and their relation to the mean number 
of photons (N) observed in a given time interval. For radiation with a degree of polarization 
< P < 1 the corresponding expression is (Wolf, 1960) 



While fluctuations can, in principle, be measured in single beams, a more reliable 
determination is often obtained by dividing the incoming radiation so half falls on one 
detector and half on another, for instance through use of a beam splitter, and determining 
the correlations or coincidence counts in the radiation incident on the two detectors. The 
correlation in the counts on the detectors, respectively labeled by subscripts 1 and 2, then 
becomes 



Hanbury Brown et al. (1967) showed that for sources small compared to a telescope's 
angular resolving power, the light could be gathered by two telescopes separated by variable 
distances. When the separation D between telescopes became comparable to the diameter 
of a telescope that, at spectral frequency v would just be able to resolve a star, the 
correlations in counts would begin to drop, giving the star's angular diameter as 6 ~ X/D, 
where A = cjv is the wavelength and c is the speed of light. While Hanbury Brown et al. 
were working with optical telescopes, following earlier observations in the radio domain, the 
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same photon bunching could also be observable with Cerenkov telescopes operating at the 
highest gamma frequencies. 

To quantitatively see how the fluctuations are related to the coincident arrivals posited 
in Paper I we need to look at the probability P(N) of finding N photons per phase cell as 
a function of mean photon number in the cell. Mandel (1963) has shown this to be 

P(A0 = (l + (N))\l + jk) N (7) 
This is the Bose-Einstein probability distribution. 

Mandel and Wolf (1995) explicitly show this expression to hold for response times r 
short compared to the correlation time of the intensity of the radiation. This is the case for 
a Cerenkov array detecting individual TeV photons arriving at intervals long compared to 
the time resolution r of the array. The meaning attached to Av may then be thought of as 
the uncertainty in the frequency of each counted photon, given by Av ~ r _1 rather than, 
as assumed by Levinson (2001), the far larger bandwidth of photons emitted by the source, 
which becomes apparent only when an entire ensemble of Cerenkov events is analyzed. 

For an aperture just large enough to resolve the solid angle subtended by the source, 
and for the fast response time of the Cerenkov array and the sparse flux of arriving TeV 
photons, tAv = 1, and g also takes on the value 1. 

Paper I has suggested that Cerenkov radiation observed at the highest TeV energies 
might arise from two or more photons arriving in coincidence from an extremely energetic 
source. It also showed how this possibility could be tested through detailed observations of 
the development of the Cerenkov shower in the atmosphere. 
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3. The Crab Pulsar and Markarian 501 

Two crucial papers now need to be cited: Aharonian et al (2000) and Aharonian et al 
(2001). 

The first of these papers reports informative tests along the lines proposed in Paper 
I. The authors find that for the Crab pulsar and the blazar Mrk 501 coincident arrival of 
photons within the angular resolution of the Cerenkov array, i.e. within a requisite width 
transverse to the direction of arrival at the altitude of initial impact, can be ruled out at 
roughly a 1 a level for the highest energy ~ 15TeV photons. The authors provide an upper 
limit on the bunching based on the anticipated development of the Cerenkov shower. At a 
90% confidence level they set an upper limit of 65% on the number of photon pairs, N = 2, 
that could be added to the shower without significantly changing the shower structure from 
that produced by single photons having an equal combined energy. At a similar level of 
confidence the authors rule out that more than 30% of the showers produced at any given 
energy could be produced by bunches of N = 10 photons, each at an energy one tenth of 
that of their combined energies. 

When Paper I was written, one could not rule out that the apparent TeV spectrum 
of a source was strongly influenced by photon bunching. The apparent spectral energy 
distribution at the highest energy could, in that case, have been entirely due to highly 
bunched photons masquerading as individual photons with much higher energies. As a 
result of the work of Aharonian et al. (2000), we now know that this possibility can be 
largely ruled out for Mrk 501, and that the observed Cerenkov spectrum well represents 
the actual spectral energy distribution of arriving photons in the ~ 1 to 10 TeV range. 
However, the concern about bunching is greatest at the highest TeV energies, where 
traversal across extragalactic space is most severely affected by annihilation through impact 
against intergalactic infrared photons. 
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The second paper of importance is that of Aharonian et al. (2001), which provides 
an improved set of spectra for Mrk 501 in both the high and the low state. Knowledge of 
the spectral shape permits us to derive an even more stringent upper limit on the photon 
bunching than given by Aharonian et al. (2000). This can be seen by considering the steep 
decline in the spectrum of Mrk 501 at the highest energies ^ lOTeV. For the year 1997, 
the observed drop in the number of photons per TeV interval was approximately a factor 
of ~ 20 in going from energies of 8 TeV to 16 TeV. If as few as 2 out of every 20 photons 
at 8 TeV arrived in coincidence, all of the 16 TeV photons could be accounted for. Since 
there are so many lower energy photons, the fraction / of these that needs to be bunched in 
order masquerade as single higher-energy photons could be quite small. A fraction / ~ 0.1 
appears quite compatible with the results of Aharonian et al (2000). From equation (3) we 
then have, (N)/g £ / ~ 0.1. 

For (N) = 0.1 the probability given by equation (7) for the arrival of an empty phase 
cell is 0.909, for the arrival of one photon it is 0.083, for 2 photons it is ~ 0.0075, and for 
each additional photon a further factor of ~ 11 less. 

A point noted by Levinson (2001) should still be mentioned. For a steady incident 
flux, Levinson set (N) / {At) equal to the flux, or the number of Cernekov counts per unit 
area per unit time, and argued that for the observed TeV flux F ~ 10 _10 cm _2 s _1 this 
corresponds to (N) £ 10~ 12 for a "spot size" 20m by 20m ~ 4 x 10 6 cm 2 and r ~ 10~ 9 s. 

This would be correct if the incident flux were steady. However, a steady flux is 
not expected. For an intermittent flux with low duty cycle, (N) can be many orders of 
magnitude higher. 

A low duty cycle is expected because the conditions that permit the production of 
gamma rays at the very highest energies are very special. Most models attempting to 
account for the TeV radiation assume inverse Compton scattering off highly relativistic 
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electrons. Paper I, for example, assumed a Lorentz factor T for a relativistic jet and a 
Lorentz factor 8 for the relativistic electrons comoving with the jet. To produce gammas at 
the very highest energies E, a product (T5) 2 ~ 10 16 (E/ (TeV), and a fortuitous alignment 
of the two directions of motion was required. The inverse Compton scatter would then 
beam these photons into solid angles of order (5^)~ 2 ~ 1CT 16 sr. The required alignment 
of the two sets of relativistic velocities, and the narrow beaming angles in what are likely 
to be turbulent jets, strongly suggest that conditions favorable to producing the highest 
energy TeV photons are short-lived and have extremely low duty cycles. The arrival of 
TeV photons will be episodic, and (N) should not be derived from a mean flux. A value of 
(N) ~ 0.1 accordingly seems well within expectations. 

Even though bunching may occur both at lower and higher energies, it is more readily 
observed at higher energies because higher energy photons from a given source are laterally 
more tightly confined within a single phase cell. If the apparent aperture, D, of the 
Cerenkov telescope — i.e. the limiting "spot size" in the upper atmosphere that can be 
resolved by the Cerenkov array — is sufficiently large to just resolve a given source at some 
energy E, then two photons at energy E belonging to the same phase cell can both give 
rise to the same shower. In contrast, as explained in connection with equation (6) above, a 
pair of photons, each having half this energy, namely E/2, will be distributed over an area 
of diameter 2D, 4 times as large as the "spot size" and the probability of both photons 
contributing to the Cerenkov pulse is reduced by a factor of 4. 

In summary, the upper limits on photon bunching that Ahraronian et al. (2000) 
provide, tell us that the observed photon number distribution N(E) is the actual spectral 
distribution, and that bunched photons masquerading as single higher-energy photons can 
only affect the spectral energy distribution to a minor extent. Given this, the rapidly 
declining photon number at high energies can be used to further constrain the number of 
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coincidentally arriving photons, simply because a small fraction of bunched photons at 
energies E/2 would suffice to provide all of the photons arriving with energies E 

Further observations, nevertheless, should be undertaken to rule out even this low 
occurrence of bunching at the highest energies because it is just these photons that are 
most critically affected by annihilation by mid-infrared photons at ~ 10 — 100 jum, for which 
direct observations of the extragalactic infrared background are particularly difficult. 

4. Fluctuations from Cosmic Lasers or Masers 

Levinson (2000, 2001) has investigated whether coincident photons could reasonably 
be expected to be produced by plane parallel lasers at TeV energies. The two papers by 
Aharonian et al. (2000, 2001) taken together now show that the observed TeV spectrum 
cannot be due to a superposition of large numbers of monochromatic photons. Accordingly, 
the flux generated at the source must be largely broad-band. This makes a laser model 
highly unlikely. 

Nevertheless, the photon statistics of laser and incoherent radiation should be examined 
in order to better understand one point emphasized in Levinson's two papers. As Levinson 
points out, the photons from a laser with diameter D must diverge with an angle ~ X/D 
because of diffraction, and the arriving beam at Earth can never be smaller than D + XL/D, 
where L is the distance to the source. For Mrk 501 at a distance of ~ 130 Mpc Levinson 
concludes that at TeV energies "the minimum spot size of a beam of photons arriving at 
Earth is of the order of a few kilometers." This is much larger than the "spot size" roughly 
20 m in diameter within which Cerenkov measurements would be able to register photons 
arriving simultaneously. 

As already mentioned, the equations presented in section 2 specifically do not apply to 
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plane parallel lasers operating well above threshold because their emission is not represented 
by a Gaussian random process, but is coherent. Under these conditions, the Bose-Einstein 
fluctuations, represented by the second term on the right in equations (5) and (6) are not 
expected. Only the first term on the right is found and takes the form of a classical Poisson 
distribution. Early experiments by Smith and Armstrong (1966) were able to verify that 
for lasers operating well above threshold the Bose-Einstein fluctuations, and with them the 
coincidentally arriving photons, disappear. 

For lasers the probability for a coincident pair of photons, as contrasted to a single 
photon, falling within the limited area viewed by existing Cerenkov arrays is then 
proportional to the ratio of the area of the array to the area of the laser "spot size" - 
roughly a factor of ~ 10~ 4 , far lower than the fraction expected for bunched photons 
arriving from an incoherent source of identical dimensions. 

The reasons for this difference are straightforward. 

The output of a plane parallel laser originates in a single spontaneous transition in 
resonance with an enclosing Fabry-Perot cavity. This seed radiation is amplified in a 
cascade of stimulated emissions that produce a beam of photons all of which are quantum 
mechanically identical. What emerges from the laser is a coherent, plane parallel beam 
effectively produced in a point source at infinity, and constricted to pass through the 
aperture of the laser end window. Such an ensemble interferes with itself to diverge in a 
diffraction limited beam that becomes wider the smaller the laser exit aperture. 

In contrast, for a source that randomly emits photons, say through blackbody emission, 
just the opposite is true. The larger the original distant source, the smaller the spatial 
region in which arriving photons with identical quantum properties are confined. This 
is because the angles of arrival, i.e. the transverse momenta of the photons are quite 
uncertain (as they would also be for a laser beam focused by a lens). This permits high 
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spatial localization limited only by the Heisenberg uncertainty principle. Processes of the 
kind advocated in Paper I, in which light is incoherently inverse-Compton scattered over a 
sizeable source area, similarly satisfy the conditions for which equations (1) to (7) apply. 

5. Discussion 

As stated in Paper I, it is difficult to know what types of astrophysical processes could 
produce bunched TeV photons in blazars, and one may argue about the plausibility of the 
specific process outlined in Paper I. Nevertheless, the very highest energy TeV gamma 
events should continue to be more precisely analyzed to determine whether a significant 
fraction is produced by bunched photons. Even if coincident arrival cannot account for 
most of the 5-10 TeV photons from Mrk 501, it will be important to place more stringent 
upper limits on photons in the 10 - 20 TeV range, where the photon flux dramatically falls 
off (Aharonian et al., 2001). If only a minor fraction of photons at energies of order 7 to 8 
TeV arrives in pairs, so that (N)/g in equations (5) or (7) is merely of order 0.1, this could 
already be sufficient to account for a significant fraction of the observed 15 TeV photons. 

Because direct observations of the extragalactic infrared background are made difficult 
by the presence of zodiacal glow, scattered starlight, and high-latitude Galactic cirrus 
emission, TeV observations may in the long run provide the most precise measure of the 
infrared background and its spectral energy distribution (Stecker, 2001; Dwek & de Jager, 
2001). For this to happen, however, we will have to be sure that the spectrum of blazars 
at different distances is well understood, and that Cerenkov measures do not slant these 
spectra toward higher energies by mistaking photon pairs of slightly lower energies for single 
photons at the very highest energies. 
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6. Conclusion 

We have shown that the two papers by Aharonian et al (2000, 2001) taken together 
permit us to restrict the number of coincidentally arriving TeV photons more precisely 
than these authors' first paper. For Mrk 501 the number of paired photons arriving at 
some of the highest energies E ~ 7 — 10 TeV cannot exceed ~ 10% without exceeding 
the total flux observed at energies 2E ~ 14 — 20 TeV, where the 1 a error bars are still 
quite large. Continued observations could, however, rule out even this eventuality. This is 
essential because the accurate determination of the mid-infrared extragalactic background 
radiation critically depends on a proper understanding of the actual energy spectrum of 
the highest-energy TeV photons. While we do not currently know of any astrophysical 
process certain to produce bunched photons at TeV energies, the possibility can only be 
observationally ruled out. 
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